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ABSTRACT 


A wind tunnel balance system is designed to determine the wind- 
induced vibrations of a space shuttle model. The balance utilizes 
a flexible sting mounting in conjunction with a geometrically scaled 
rigid model. Bending and torsional displacements are determined 
through strain- gauge- instrumented spring bar mechanisms. The natural 
frequency of the sting-model system can be varied continuously 
throughout the expected scaled frequency range of the shuttle vehicle 
while a test is in progress by the use of moveable riders on the 
spring bar mechanism. Through the use of a frequency analyzer the 
output can be used to determine troublesome vibrational frequencies. 

A dimensional analysis of the wind-induced vibration problem is 
also presented which suggests a test procedure. In addition a computer 
program for analytical studies of the forced vibration problem is 
presented . 
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FOREWORD 


This document presents the results of work performed at the 
Louisiana State University, Baton Rouge, under the Graduate 
Engineering - Practice in Mechanical Engineering Contract, Number 
NAS9- 10464. The report covers the work performed during the 
Spring and Summer semesters, 1971. 

During the Spring semester full-time student participants were 
Mr. T. Doran and Mr. R. Eslami, while part-time participants were 
Mr. G. Plaisance and Mr. J. Rubli. The latter individuals performed 
preliminary work on the project while completing their own work in 
the program on other topics which was begun at an earlier time. 
Messrs. Doran and Eslami continued on the present project through 
the Summer semester. Throughout the reporting period the Program 
Director was Dr. L. R. Daniel, the Associate Director was Dr. R, W. 
Courter, the Principal Advisor was Dr, P. H. Miller and the NASA/MSC 
Project Monitor was Mr. Charles Teixeira. In addition, during the 
Spring semester a graduate course on random vibration of structures 
was taught by Dr. G. D. Whitehouse. 
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INTRODUCTION 


The development of a space shuttle system is currently in 
progress under the guidance of the National Aeronautics and Space 
Administration. This program involves not only the development of 
new technology, but the extension of existing technology to broader 
scope. The problem of wind- induced vibrations on shuttle vehicles in 
the launch configuration is expected to be more severe than it has 
been on launch vehicles used to date. Because of the complexity of 
the shuttle structural arrangement and the non uniformity of the 
aerodynamic field, the determination of critical vibration conditions 
may prove to be a difficult and certainly tedious task. 

In order to simplify the design process it is advantageous to 
devise a method whereby critical areas of wind- induced vibrations 
can be immediately and inexpensively determined. This report is 
concerned with the analysis and design of a wind tunnel balance 
system which will perform this task. 

The report which follows is divided into two parts. Part I 
deals with the specific details of balance system design and Part 
II describes some analytical work which was performed to determine 
instrumentation requirements and to develop an analytical capability 
for the study of experinental results to be provided by the balance 
during a test program. 
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CHAPTER I 


INTRODUCTION 

When vehicles are standing on the launch pad prior to launch, 
they are often subjected to dynamic loads from prevailing winds and 
gusts. In the past, some of these dynamic loads have been suffi- 
ciently severe to produce dangerous vibrations In the launch vehicle. 
In these cases. It was necessary to apply vibration dampers to the 
vehicle at the launch site to prevent dangerous vibration levels. 

For years, engineers have grappled with the problem of wind 

loads on flexible bluff bodies. Many problems develop In the 

determination of wind loads on such structures as smoke-stacks, 

'ic 

tall masts and suspension bridges [1] . A similar problem of the 
determination of b&th static and dynamic wind loads Is encountered 
with the shuttle. 

In recent years, much study has been done on vibrational 
characteristics of launch vehicles. Much study was done on either 
full or large scale models, and with mathematical models (finite - 
element program NASTRAN) [2]. Tests of this sort require either a 
"frozen" design or the actual vehicle on the launch pad. This Is 
not possible with the shuttle since design Is not complete and, 
of course, no vehicle has been constructed. No test system now In 
use Is applicable to both fueled and unfueled shuttle. The unfueled 
shuttle Is considered to be more sensitive to wind Induced vibra- 
tions than Is the fueled vehicle. It Is desired that not only 

* 

Numbers In brackets Indicate references listed on page 
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steady-state wind loading but also gustlng and variable wind profiles 
be studied [l]. 

It Is expected that ground wind loads will be even more 
critical for the unsymmetrlcal space-shuttle configuration than 
they have been on near- symmetrical vehicles which have already been 
launched. This is true because of the non-uniform, non-circular 
cross-section, large lifting and control surfaces of the shuttle 
and the expected "piggy-back" arrangement of the orbiter on the 
booster. 

Several possible types of induced motion can be identified: 

1. "stop-sign" vibration 

2. "galloping" vibration 

3. transverse bending vibration 

"Stop- sign" (torsional) vibrations are caused by such phenomena 
as wing stall. "Galloping" (longitudinal bending) vibrations are 
caused by vortex shedding at different frequencies along the body 
of the shuttle because of variations in cross-section. Transverse 
bending vibrations are caused by the shedding of vortices in different 
directions along the body. The input which creates these vibrations 
is a function of wind velocity, angle of wind direction with respect 
to shuttle, and configuration of vehicle. Other factors which affect 
the aerodynamic input are wind-profile, ground plane shape and 
proximity of service towers which may create turbulence in wind. 

It is desired to design a test system which can be used in a 
wind tunnel to evaluate the effect of ground winds on the launch 



pad configuration of the shuttle. The test system is to be a small, 
inexpensive system which would be applicable to various shuttle 
configurations . 

The simulation requires selection of important simulation 
parameters and modeling of the torsional and bending modes of 
vibration. This is accomplished by simulating the flow field with 
a geometrically scaled rigid model. The fundamental response of 
the model is simulated through accurate inertia scaling and mounting 
on a flexible sting which will produce scaled bending and torsional 
vibration frequencies. It would also be possible to simulate 
ground plane wind profile, and service towers and other physical 
structures which may influence flow. 

This test system will be used as a first test to define areas 
of critical vibrational loads. Later tests will use larger models, 
probably not rigid, to further study these critical conditions. 
Considerable time and money can be saved by this test system since 
much less time and expense will be needed for the study of non- 
critical areas with more elaborate test equipment. 

At the present stage of design the shuttle payload is 
approximately 1% of the total vehicle gross liftoff weight. If 
areas of critical conditions can be related to structural design 
early in the vehicle design, it could be possible to produce a 
weight saving while avoiding critical vibrational conditions. If 
wind tunnel testing of the shuttle in ground winds is delayed until 



the shuttle design is complete, the correction of discovered initial 
conditions could require a weight penalty which could actually 
estimate the entire payload [3]. This indicates just how weight 
critical is the entire vehicle. 



CHA.PTER II 


SCALING CONCEPTS 


In all wind tunnel modeling It is necessary to equate certain 
dimensionless quantLtlfis to insure that data from the model in applic- 
able to the prototyp>e. In tests where wind-induced vibration is 
involved, as in flutter experiments, the frequency of input must be 
scaled by a dimensionless parameter. In modeling the shuttle in 
launch-pad configuration, it is necessary to model both the flow, 
and the vibrational characteristics. 

Since it is desired to determine critical vibration conditions 
on the full scale vehicle by testing a rigid model mounted on a 
balance in the Texas A 6e M Low Speed Wind Tunnel facility, the model 
must physically fit into the test section of this wind tunnel. The 
test section is 7 feet by 10 feet wide by 16 feet long and operates 
at room temperature. Maximum free-stream velocity in the test sec- 
tion is 300 ft. /sec (204.5 mi/hr). [4] It was concluded that a 
scale factor of 1/70 would create negligible blockage. This is 
considered smaller than the largest scale factor that could be used 
without introducing large error due to blockage and wall effects. It 
is desired to use a scale factor as large as possible to allow greater 
accuracy in model dimensional detail. Scale factor is defined as: 


SF 


Dp 70 


(II-l) 


5 
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Following selection of scale factor, other dimensionless para- 
meters to Insure similitude between model and prototype are to be 
found. 

In most wind tunnel tests, it is desirable to have similar 
Reynolds Number, defined as: 


Re “ 


p V D 


(n-2) 


If Re 3 molds Number for model and prototype are equated, then 


Re^ “ Rep 
or 


(II-3) 


D V D 
M M 


M 


Pp '^p Op 


(II-4) 


Since air is to be the common medium: 






(II-5) 


It is noted that either compressed air or freon gas may be used in 
some wind tunnels to facilitate simulation of Re, but the Texas A 6e M 
wind tunnel facility has no capability for using compressed air or 
freon gas medium. 

Equation (II-4) reduces to 






(n-6) 
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Solving for 


yj a# V ■“ 

M p n 


M 


Substitution of Equation (II-l) into (II-7); 


(II-7) 


M 


P SF 


(70) 


(II-8) 


Since the maximum wind velocity to be experienced by the shuttle 

at Cape Kennedy is 70 mi/hr, this will be the maximum value of Vp. 

[ 5 ] Substitution for v„ in Equation (II-7) yields: 

P max 


Vj^ = (70 mi/hr). (70) = 4900 mi/hr (II-9) 

This result leads to the conclusion that a Reynolds Number 
similitude is not possible in a low speed wind tunnel. Low Reynolds 
Number tests have been done by NA,SA/Langley and research is now under- 
way to correlate these data to high Reynolds Number condition. [ 2 ] 

The next dimensionless parameter to be investigated is Mass 
Ratio, defined as: 


MR = 


P ^ D S 
^air 

m 


If Mass Ratio is to be similar, then: 


(II-IO) 


MR^ » MRp 

Expanded this becomes: 


(II-ll) 


air 


°M ^ 


"m 


p D S 
^air P P 

m„ 


(n-12) 
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Simplification gives; 

m D S 
-jE » ^ _P 

“m ^ 

A case of Equation (II- 1) says: 


(H-13) 




S 


P 


(SF)^ 


(11-14) 


Substitution of Equations (II-l) and (11-14) into Equation (11-13) 
yields : 


“m (SF)^ 


( 70 )^ 


It is known that 


m = p V 


(n-15) 


(11-16) 


Substitution of Equation (11-16) into (11-15) yields: 


M 

P 



(70)^ 


Again a case of Equation (II-l) gives: 


(n-17) 


O -I 


(11-18) 
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Therefore , 



Pp 


(SF)^ 


(SF)^ “ 1 


or 


P 


M 


P 


P 


(n-19) 


(n-20) 


Equation (11-20) essentially states that for similitude of Mass 
Ratio, the average density of the model must equal the average den- 
sity of the prototype. 

Also of Interest in dynamic modeling is Inertia Distribution, 
Inertia Distribution is defined as: 


ID 



Similitude yields: 


This may be written: 



5 5 

Solving for yields: 

5 



(n- 21 ) 


(11-22) 


(11-23) 


(11-24) 
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Since for similitude of Mass Ratio: 


P 

P 


M “ 1 
P 


(11-19) 


Then 





(11-25) 


Equation (11-25) can be reduced by use of Equation (11-21); 


5 

= Jp . (sr)= - Jp (ij) (11-26) 

This equation is used to model mass moment of inertia. 

Probably the most important and widely used parameter for vibra- 
tion modeling is Strouhal number. Strouhal number is widely used in 
flutter experiments to relate input frequency to wind velocity. [?] 

It is defined by; 


St= (H-27) 

Although frequency in reference to Strouhal number refers to 
aerodynamic input frequency, in this design case, it will be used 
to represent the natural frequency of the model or prototype. This 
assumption is justified by the fact that the frequency of the aero- 
dynamic excitation equal to the natural frequency of the shuttle is 
the critical conditions to be studied. 

Similar Strouhal number between model and prototype yields: 
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S*=M “ 

or 


(11-28) 


\ \ \ 

Solving for : 

M 

f = f i 
\ ^ °M "p 

Substitute Equation (II -1) gives; 




(70) 



(11-29) 


(11-30) 


(11-31) 


It is noted that in Equation (11-31), the ratio ('fj^/Vp) is vari- 
able, that is, it has not been specified by any other dimensionless 
parameter similitude. Of course, by wind tunnel specification, is 
limited to a maximum of 204 mi/hr. As mentioned earlier, the maximum 
wind velocity impinging on the shuttle is 70 mi/hr. From the maximum 
wind tunnel velocity and the maximum Cape Kennedy velocity, it is 
seen that the velocity ratio defined as: 

V 

VR = — (11-32) 

P 

has the range 


0. < VR < 3. 


(11-33) 
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It is noted that a higher Velocity Ratio will increase output 
amplitude by increasing static wind loads. However, this has the 
effect of raising model natural frequency, which as will be explained 
later is not desirable. Therefore, a Velocity Ratio is chosen: 

VR = ~ = 1 (11-34) 


Substitution of Equation (11-34) into Equation (11-31) gives: 


Equation (11-35) is used to model natural frequency of the test 
system. Data obtained from NA.SA./Langley indicates a range of natural 
frequencies for the shuttle of 0.2 to 4.0 Hz will encompass the first 
four modes of vibration for the combined fueled booster and orbiter 
in launch pad configuration. A range of natural frequencies from 
0.8 to 8.0 Hz will encompass the first four modes of vibration for 
the unfueled orbiter and booster combination. [6] 

It is possible to model the shuttle in launch configuration by 
use of the following modeling parameters and respective modeling 
equations: 

1. Scale factor 



1 _ 

70 


(II-l) 
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2. Mass Ratio 


P = p (avg)p 

3. Inertia Distribution 


(11-20) 


= Jp (SF)- 


( 70 )^ 


4. Velocity Ratio 


(11-26) 



5 . Strouhal Number 


(11-34) 





. (70) 


(H-35) 


These parameters are to be used in design of the test system. 



CMPTER III 


ESTIMATION OF WIND LOADS 

After selection of similitude parameters is complete, it is 
necessary to be able to estimate maximum wind loads. The magnitude 
of wind load does not affect scale and sizing of model, which are 
set by similitude parameters, but does affect the balance members, 
the construction and the materials incorporated in the design of the 
balance. 

Aerodynamic input is very difficult to predict, even on very 
simple geometric shapes such as cylinders. The frequency and ampli- 
tude of the input from vortex shedding. [3] 

Being random, the aerodynamic input signal and its maximum 
amplitude and impossible to describe analytically. This is precisely 
why a wind tunnel test is needed; that is, to determine the frequency 
content and other characteristics of this random aerodynamic input. 

It is concluded that the static wind load would represent a 
large part of the input, and could be roughly calculated with 
known data for similar configurations. 

It is possible to determine the forces on a body if the 
wind velocity, coefficient of lift and drag, density of impinging 
medium, area and area centroid of the body are known. Any body can 
be broken down into a finite number of simple geometric shapes and 
a summation taken over all parts to approximate the total reaction 
on the body. 
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The shuttle configuration now under study by NA.SA/MSC is a 

delta-wing orbiter mounted "piggy-back” on a delta-wing booster 

2 

(Figure III-l). The booster has a planform area of 13,439 ft, 

2 

and the orbiter has a planform area of 8,442 ft. . [8] Planform 

area of the combined vehicle was not available but was graphically 

2 

estimated to be 20,000 ft. . 

The planform view of the shuttle was divided into six sections 
with simple geometric shapes (Figure III-2) with easily found centroids 
(centers of pressure) of each of these sections. Dynamic pressure 
is defined as [9]: 

Q = % (III-I) 

and the lift and drag forces on each section is 

Drag = Cjj . A . Q (IH-2) 

Lift = . A . Q (IH-3) 

The components of lift and drag are given in Figure III -3 as 
related to planform model. These are not exact since the mutual inter- 
ference is not taken into account. 

As the vehicle will be subjected to winds from all directions, 
it is necessary to perform the static wind load calculation over a 
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G COMPONENTS 


III-3 
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range of 180° (only angles of attack of 0 ° to 180° must be studied 
since the shuttle is S3msnetric) . 

It is noted that in Equations (III- 2) and (III-3), and 
are functions of angle of attack. However, values of Cj^ and 
for two configurations, fin- tip orbiter on fin- tip booster, and 
straight-wing orbiter on delta-wing booster (Figure III-4) , are 
known as a function of angle of attack. [3] It is considered that 
the values of C and C for the delta-delta shuttle would represent 
no drastic change in C or C from either of the above configurations. 

Using the above values for and and planform model in 

Figure III-2, it is now necessary to perform the calculation of the 
three moment defined in Figure (III-5). The equations of those 
moments are: 


T = (D3 cos 0 + L3 sin 0) • (SL4/2 + SL3/3) 

+ (Dl cos 0 + LI sin 9 ) • SL4/4 

(III-4) 

+ (D2 sin (0-cp) + L2 cos (0-cp) • 5*SL4/6 sin cp 
+ (D4 sin 9 + D4 cos 0) • SL/1/2 

ST 1 

LB = (Dl sin 0 - LI cos 0) • + TC 

CT 2 

+ 2(D3 sin 0 - L3 cos 0) • ^ + TC 

4. ST 2 

+ 2(D2 sin 0 - L3 cos 0 ) . + TC 

ST 1 

+ (D4 Sin 0 - L4 cos 0) • + TC 


(III-5) 





DEFINITION OF MOMENTS 


figure IH-5 
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TB = (D1 cos 0 - LI sin 0) 
+ 2(D3 cos 0 - L3 sin 0) 
+ 2(D2 cos 0 - L2 sin 0 ) 
+ (D4 cos 0 - L4 sin 0 ) 


2 

SL2 

3 

4SL2 

3 

SLl 

2 


(III- 6) 


where SLl = 284 Ft * Scale Factor 
SL2 = 90 Ft * Scale Factor 
SL3 = 60 Ft * Scale Factor 
SL4 = 30 Ft * Scale Factor 

as illustrated in Figure (III-3). A program (Appendix A) was written 
to perform these calculations for each 7.5° increment of angle of 
attack. Results from this program are listed in Appendix A. The 
maximum expected static moments are for fin tip shuttle: 


T =8.35 ft- lb 

(III- 7) 

max 

LB = 81.35 ft- lb 

(III- 8) 

max 

TB = 96.92 ft- lb 

max 

(III-9) 


And for the straight-delta configuration: 


T =8.50 ft- lb (III- 10) 

max 

LB = 84.65 ft- lb (III- 11) 

max ' 

TB = 92.72 ft- lb (III- 12) 

max 

The results of the program are presented graphically in Figures 
III- 6, III- 7, III- 8. 



fin-tip configuration 



MAGNITUDE OF TORSIONAL MOMENT ANGLE OF ATTACK (Deg.) 

FIGURE III- 6 




MAGNITUDE OF LATERAL BENDING MOMENT ANGLE OF ATTACK (Deg.) 

FIGURE III- 7 



straight, delta-wing configuration 



MAGNITUDE OF TANGENTIAL BENDING MOMENT ANGLE OF ATTACK (Deg. 
FIGURE III-8 



CHAPTER IV 


PRELIMINARY DESIGN 

Before final design of the test system is attempted, it is 
necessary to set down certain properties which the test system 
must possess. Most important, the test system must insure 
simulation of all similitude parameters discussed in Chapter II. 

The test system is divided into three sections: 1) the wind 

tunnel, 2) the rigid model of the shuttle, and 3) the flexible 
sting which will simulate vibrational properties of the shuttle. 
The wind tunnel will be used to simulate velocity ratio 

''m = ''p (”-34) 

The model itself will simulate mass ratio 

^avgM " ^avgp 

From Ref. [7], the total mold line volume is given as 
= 25,671 ft^ 

c 

Vq = 123,168 ft^ 

and the mass of each fueled vehicle is 

nig = 3,936,387 Ibm 

m^ = 944,804 Ibm 
F 


(II- 20) 

(IV- 1) 
(IV- 2) 

(IV- 3) 
(IV-4) 
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and unfueled: 


= 633019 Ibm 

®UF 

(IV- 5) 

m^, = 276115 Ibm 

(IV- 6) 

UF 


average density, P^^g* found by: 


- E 
^avg V 

(IV- 7) 


Then it may be concluded that 


p = 15.33 Ibm/ft^ 
“f 

(IV- 8) 

p = 2.47 Ibm/ft^ 

%F 

(IV- 9) 

p =7.67 Ibm/ft^ 
F 

(IV- 10) 

p =2.24 Ibrn/ft^ 

UF 

(IV- 11) 


Models of the above four densities can be made either of 

3 

balsa wood (density 7.49 - 12.79 Ibm/ft ) or polyurethene foam 
[10]. These models must be geometrically scaled to the scale factor 



(II- 1) 


For scaling of inertia distribution it is necessary to build 
model and balance such that 


■^TOTAl" ‘^MODEL '^BALANCE “ 


(IV- 12) 
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From Reference [8] the mass moment of inertia for the axis 
system in Figure III-5 is given as; 

Jp = 33.507 X 10^ slug-ft^ (IV- 13) 

XX 

Jp = 530.434 X 10^ slug-ft^ (IV- 14) 

yy 

J = 515.492 X 10^ slug-ft^ (IV-15) 

^zz 

for the fueled vehicle about the center of gravity. 

Substituting into Equation IV- 12 and solving for J about base 
of booster (point 0, Figure III-5) gives: 

Jtotal =0*6^2 Ibm-ft^ (IV- 16) 

XX 

Jtotal “ Ibm-ft^ (IV-17) 

yy 

Jtotal = ^o. 12 ibm-ft^ (iv-i8) 

zz 

Reference [8] does not give J for unfueled booster and 
arbiter, but does give J for total vehicle at booster burnout. It 
is concluded that these data are not very different from those 
for the totally unfueled vehicle. Therefore these values were 
used to calculate "ball-park" figures for J of the balance for 
unfueled configuration. From Reference [7]; 

Jp = 28.781 X 10^ slug-ft^ (IV- 19) 

XX 

Jp = 229.870 X 10^ slug-ft^ (IV- 20) 

yy 

Jp = 221.194 X 10^ slug-ft^ 

zz 


(IV- 21) 



for the vehicle at booster burnout at the center of gravity. 
Relating these data to equation IV- 12 and moving the axis system 
to point 0 (Figure III-5) gives: 


^TOTAL “0*55 Iboj-ft 

(IV- 22) 

XX 

'total = 

(IV- 23) 

yy 

TOTAL - ”-35 Ibn^ft^ 
zz 

(IV- 24) 


for the unfueled vehicle model. 

The balance itself will be used to simulate Strouhal number: 

• (70) (II- 35) 

The balance is to be designed so that the model (mounted on 
the sting) has 3 degrees of freedom. The three motions (similar 
to Figure III-5): are l)Torsion, 2) Lateral Bending, and 3) 

Tangential Bending. All three motions are about the point 0 
(Figure III-5). 

The equation of a system oscillating about a point is [ll]: 

+ C0 + Ke = F(t) (IV- 25) 

where 0 = position angle 
C = damping 

K = spring constant ft- lb /rad 


F(t) = forcing function 
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If damping Is neglected: 


+ K0 » F(t) 


•• 2 F 

9 + (Ujj 0 = j (t) 


where 


K ,rad , 
J '‘sec' 


Here, is the natural frequency of the system. We know that 


ou = 2n f 


Combination of equations (III- 35), (IV- 28) and (IV- 29) gives 


From Ref. [6] ranges of natural frequencies of the fueled and 
unfueled vehicle are given: 


0.2 Hz ^ f„ ^ 4.0 Hz 


0.8 Hz ^ f„ i 8.0 Hz 
^UF 


Combining equations (IV- 28) and (IV-29) gives: 




(IV- 33) 



Solving for K; 


K = (2TTfj^)^ • J (IV- 34) 

From this equation maximum and minimum values of K can be 
calculated. Different values of K must be calculated for torsion 
and bending, since J for each of these motions is different. 

For torsion: 


K, = 1700 ft- lb /rad 
low 

= 1,104,500 . ft-lb/rad 
and for bending 


^low 

^high 


= 202131. ft-lb/rad 
= 115,811,500. ft-lb/rad 


It is desired for the balance to have the capability of 
variable frequency in the specified range. This property would 
facilitate use with any natural frequency and would allow testing 
of a configuration for any natural frequency. With a variable 
frequency balance, tests to determine frequencies to be avoided 
could be easily done. 

Since the torsional and bending vibration may be coupled, 
it is desired to design a balance which will allow testing 
simultaneously. 



The desired output of the test system is an energy vs. 
frequency graph for a given natural frequency of the vehicle. It 
would be advantageous to have a direct method of instrumenting 
this graph directly from the test system. 



CHAPTER V 


DETAIL DESIGN 

The most Important property of the balance is a variable fre- 
quency capability. For a constant mass moment, J, this is accomp- 
lished by varying the spring constant, K, of the system. In this 
design this is accomplished by changing the length of a flexible 
bar ("spring bar"). 

The spring constant for this type of spring must be derived. 
For a bar with moment M (Figure V-1), the slope at A is: [12] 


e 


1 ^ 

3 El 


(V-1) 


In general, the equation for a torsional spring is: 


M = K 0 (V-2) 

Rearranging Equation (V-1) gives: 

M = 0 (V-3) 

Jj 


frcm which can be seen that 

1j 


(V-4) 


for this type of loading. 
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The design Incorporates two such "spring bars" in parallel, 
therefore : 

■Sotal - (''-=> 

For a given material (E) and a given area moment (I), the 
spring constant can be changed by changing the length of the bar. 

This is the method to vary K in this balance. 

A program was written to size the "spring bars" for steel 
(E = 30,000,000 psi) and a circular cross section (I = tt R^/4). 

This program and results are listed in Appendix B. 

This type of "spring bar" is used to model both torsional and 
bending modes. Three separate bars are required to cover the entire 
bending frequency range desired. Five different bars are needed to 
cover entire range of torsional frequencies. 

The torsional balance is made with five "spring bars" to facilit- 
ate reducing its size so that it can be built into the bending 
balance without adding too large a mass moment to the balance. 

It is desired to keep the natural frequency of the balance as high 
as possible to prevent second mode effects of balance and model. 

Complete views of balance are shown in Figures (V-2) to 
(V-21). 

It is necessary to determine the mass moment for the balance 
about an axis system passing through the universal joint. To do 
this, the balance was divided into 16 sections with simple geo- 
metric shapes and the mass moment about the respective axis system 
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COMPLETE BALANCE SYSTEM 
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TORSIONAL BALANCE 
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BENDING BALANCE 
FIGURE V-4 



BALANCE (SIDE VIEW) 





-STING 
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BALANCE (SIDE VIEW) 





"to® 
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TORSIONAL BALANCE (TOP VIEW) 
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BALANCE CENTER BLOCK (TOP VIEW) 
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BALANCE UPPER SECTION (TOP VIEW) 
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DIMENSIONS FOR MOVEABLE RIDER 


TABLE V-1 
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DIMENSIONS FOR MOVABLE RIDER 
(TORSION) 


D 

R 

1 

.1325 

I 

. 1575 

15 

.1925 

15 

.1975 

20 

.2425 

30 

.2825 

30 

.3125 

35 

.3675 


TABLE V-2 




calculated for each section. A program (Appendix A) performed these 

calculations. The results are: 

\ 

J = 0.0243 16 in/ft^ 

XX 

J = 1.37 16 in/ft^ 

yy 

J = 1.27 16 in/ft^ 
zz 

From this it is determined that the mass moment of the model 
above must be 

J = 0.640 - 0.024 = 0.616 1$ in.ft^ 

XX 

J = 70.41 - 1.37 = 69.04 16 in.ft^ 

yy 

J = 70.12 - 1.27 = 68.85 16 in.ft^ 
zz 

It is recommended that balance be tested after it is built to 
obtain the exact value of mass moment for each axis or a "shaker" 
can be attached to the model to determine its response to a sinu- 
sodial input. 

The balance is "tuned" to a natural frequency by moving the 
riders either in or out, or changing the "spring bar". The riders 
are moved and their position known by the screw drive and the 
number of turns the screw drive bar made. 



The rider slide in dovetail track and contact the spring bars 
with rollers, top and bottom. The roller permits essentially no 
axial forces to exist in the spring bars. 



CHAPTER VI 


INSTRUMENTATION AND DATA - REDUCTION 

In any wind tunnel test, instrumentation and data reduction are 
very important. It is necessary for a successful test procedure to 
be able to monitor certain outputs and apply these outputs to 
determine the desired output of the tests. If necessary, it is 
possible to record results and analyze them on a digital computer or 
other device. 

Since this test system is to be very simple, fast and inexpen- 
sive, use of the digital computer for data reduction may be avoided. 
The most desirable system is an instrument package which will 
directly output desired information. 

The requested output by NASA/MSC is an energy vs. frequency 
curve for a given natural frequency and configuration. It is 
possible to obtain this output directly in this test system. 

The instrument package consists of four strain gages, mounted 
on each moveable rider. Temperature compensating gages are mounted 
to base plate. A voltage divider circuit is set up using the gages 
in series as one resistance. A constant voltage source is required. 

As a load is put on the rider, the resistance of the gages 
changes and the voltage across the gages increases and this change 
in voltage is amplified by a DC amplifier (frequency response: near 

DC to 1 K Hz). The output of the amplifier is the input to a wave 
analyzer (frequency dependent VTVM) and this output is recorded on 
a strip-chart recorder. The graph recorder on the strip chart will 
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be a power spectral density (PSD) of the vibrational (or response) 

of the vehicle. The PSD will relate the amount of power (or energy) 

at each frequency for a given natural frequency. The units of 
2 

PSD are x /Hz as plotted against 

Also to be plotted is PSD of the acceleration of the model. 

This is accomplished by mounting piezoelectric accelerometers in 
the model. Output of accelerometer is connected to a change 
amplifier and from the amplifier to a wave analyzer with strip-chart 
recorder. 

Acceleration and displacement are to be measured and recorded 
in the three directions of (Figure III-5), i.e., torsion, lateral 
bending, tangential bending. For each run, six graphs will be made. 

Runs must be made for each wind velocity angle , 0 , wind 
velocity, v, and natural frequency, f^^. If the angle, 0, is 
incremented by 7.5° and velocity is incremented by 5 mi/hr and 
the frequency is incremented in 0.2 Hz intervals, approximately 
7000 runs will be needed to completely analyze one configuration. 
Many of these runs may be deleted for various reasons, such as for 
velocities below 30 mi/hr. This would cut total number of runs to 
5000. Figuring 10 minutes per run, total test time is approximately 


700 hours. 



CHAPTER VII 


CONCLUSIONS 

A wind tunnel test system has been designed that will facilitate 
testing the effects of ground winds on the space shuttle vehicle in 
the launch configuration. The test system includes a balance system 
and an instrumentation package. 

The balance system is used to model the fundamental vibrational 
response of the shuttle and is applicable to various shuttle config- 
urations. The balance will simulate the first four modes of vibration 
of the shuttle in either the fueled or unfueled configuration. The 
balance was designed to allow testing of the torsional and bending 
modes simultaneously since the vibrations may be coupled. 

The instrumentation package has been designed to provide a direct 
output of the Power Spectral Density of displacement and acceleration 
for the response of the vehicle to wind induced vibration. This 
closely approximates the energy vs frequency plot desired by 
NASA/Langley and NASA/MSC. 
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IZ « MASS moment around Z AXIS (LBM>FT«*2). 

IXX = TOTAL MASS MOMENT AROUND X AXIS (LBM-FT**2) 
IZZ = TOTAL MASS MOMENT AROUND Z AXIS <LBM-FTA*2) 
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MOMENT = 0.64 LBM-FT«*2 OIA.fSPRlNG BAR) = 0.140 IN. MAX TORQUE 120. IN-LB 
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2.6 14.963 0.214 0.002436 2B.57 
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MASS MOMENT = 0«64 LBM-FT**2 01 A.( SPRING BAR) = 0.190 IN. MAX TORQUE = 120. IN-LB 
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4*2 300*320 4*290 0*008676 203*39 

4*3 291*859 4*169 0*008393 196*72 

4*4 284*074 4*058 0*008127 190*48 

4*5 276*882 3*955 0*007877 164*62 
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PART II 


ANALYTICAL STUDIES 
by 


Reza Eslami 



An Analytical Method for Wind- Induced Vibration Prediction 


A dimensional analysis was made to predict and model the natural 
frequency in bending and torsion. It was found that the natural 
frequencies depend upon only the structural stiffnesses and the 
moment of inertias of the system in bending and torsion. These 
relations are as follows 


n^^(x) = C^(x) 


K^.(x) 


9 K^(x) 

■ nb <x) = ^00 

b 


where 


X = Coordinate axis along the vehicle 
n^(x) = natural frequency in torsion 
Cj^(x) = some constant 
K^(x) = torsional stiffness 
I^(x) = torsional moment of inertia 


n^,(x) 

C2(x) 

Kb(x) 

IbU) 


natural frequency in bending 
some constant 
bending stiffness 
bending moment of inertia 


( 1 ) 

( 2 ) 


[ L ] 

[T-^] 

[ 1 ] 

[ML^] 

[t"^] 

[ 1 1 

9 -9 

[ML T ] 
[ML^] 
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Equations (1) and (2) are written in terms of coordinate x. 

Thus they may be represented by continuous or step-wise functions of 
X. In the present design analysis the whole system has been 
lumped into one step and therefore one degree of freedom. Thus, 
the equations (1) and (2) could be simply 


n 


t 


2 



2 

"b 



(3) 

(4) 


in which resemble the characteristics of the whole system. Writing 
equations (3) and (4) for the model and the prototype give 





(5) 


( 6 ) 


where primes refer to the model and unprimes to the prototype. 

With the assumption that the system is just one degree of 
freedom the analytical solution of the motion of the system is 
possible. The governing differential equations of the system are 
as below 


T •• 

It'P 


+ c. 


9 


+ K cp = M 


aero 


(7) 


^b^ + ^ - ""a 


( 8 ) 
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where cp and y are torsional and lateral displacements, respectively, 

as shown in Figure 1. In the absence of the gantry and mechanical 

damper, terms and in equations (7) and (8) will drop out. The 

forcing functions, the equations (7) and (8) 

are assumed to be constant and are obtained^. With these assumptions 

the solution of the equations (7) and (8) will become 

X 



= moment of inertia of the model 
and balance system around 
the point 0 


= moment of inertia of 
the model and balance 
system along the axis 
x-x. 


K. => bending stiffness of the bars AC 
and CB 


= stiffness of the 
torsional sensitive 
element 


1 


Progress Report of May 3, 1971. 
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M 


cp(t) = sin X^t + C 2 cos X^. 


t 4* 


aer o 
L + D 


y(t) = sin X^^t + cos Xj^t — 


a 


where 



(9) 

( 10 ) 


( 11 ) 


where C^, and are the constants of the integrations. 
Assuming a typical initial condition as follows 


t = 0 t = 0 

cp = cp^ rad. (12) 

y = y^ cp=y = 0 


Substituting the initial conditions (12) into the above equations 
yield the displacement of the center of gravity of the system as 

M M 


cp(t) = (cp^ 


aero^ , . . 

— ) cos X. t + 

L + D 


aero 

^t 

L + D 


(13) 


y(t) = (y + — JT"^) 


cos x^^t - 



(14) 


The above equations will simply show that motion of the model 
is a plane harmonic function in which its frequency is the natural 
frequency of the system. 
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An investigation was made on using the Sandia shock Code to 
predict the motion of the model, but since this computer program 
will just solve lateral and axial problems the work was stopped. 

At this stage, an attempt was made to write another separate 
computer program which would more-or-less do the same calculations 
as shock code, but for lateral and torsional vibrations. The resulting 
computer program basic ly uses the fourth order Runge-Kutta method 
and predicts the displacement, velocity and acceleration of the 
system. This program represents a capability for solving coupled 
lateral and torsional vibration problems with forcing functions 
of any form (including random). It can be used to investigate 
shuttle response characteristics for various aerodynamic input 
functions if later investigations of this nature are required. 

For the present work response to simple functions of the one 
dimensional system can be adequately treated analytically without 
resorting to the numerical solution. This work is performed here 
to investigate the instrumentation requirements for the proposed 
balance system and to get a sample response for a periodic input. 

As before, the corresponding differential equation of the system 
can be represented by 




+ C cp + K 


C^cp = 


aero 


(15) 


T- •• 

^b^ 






L - D 

a 


(16) 
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which are essentially equations (7) and (8). Again it was assumed 
that the damping term is zero. The aerodynamics moment and lift 
and drag are assumed to vary according to a sine function with an 
angular velocity of u). This is similar to the Von Korman aero- 
dynamic forces. The intensity of this alternating force can be 
written as: 

p f A (17) 


2 

where is the Von Karman force, p V /2 the stagnation pressure, 

A the projected area of the obstacle (perpendicular to the stream), 
and a dimensionless coefficient, which may be called the Karman 
coefficient. For a cylindrical body the maximum intensity of the 
aerodynamic forces (such as lift, drag or aerodynamic moment) can 
be written as 




A) s in o)t 


(18) 


The same form of force can be assumed for the lift, drag and 
aerodynamic moment. Thus 


M = M sin cut (19) 

aero o 

L + D = (L + D) sin cut (20) 

where (L + D) is the resultant of the lift and drag. Substituting 
into the equation of motion gives 
_ 


Reference by Den Hartog is attached to the report. 



I^cp + K^cp = M sin cot 
Ij^y + K^y = (L + D) sin cot 


( 21 ) 
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( 22 ) 

Assuming that the initial conditions are as follow 


t = 0 

t = 0 


II 

o 

II 

-s • 

II 

o 

(23) 

9 = 'Po 




yields the following solutions for the motion of the center of 
gravity of the system 

M /K 

cp(t) = cp cosX t 5 — — (-^^ sin X t - sin cot) (24) 

° (x/-co^) ^t ^ 

(L+D )/K 

y(t) = y„ cosx. t 5 — - — sin X, t - sin cot) (25) 

o b > 2 2,. X. b 

(Xfo -u) ) D 

The velocities are then 

coMo/Kt 

cp(t) = - cp^Xj. sin Xj.t 2 — | (26) 

Xj. -CO 

co(L+D )/lC 

y(t) = - y^X|^ sin X^^t — ^(cos X^^t - cos cot) (27) 

Xj^ - CO 

and finally, the accelerations are 

$(t) = (M^ sin cot - K^.cp)I^. (28) 


y(t) = [(L+D) sin cot - K^y]/I^ 


(29) 
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where, in the above equations and X^ are the same as equations (11). 

From the equations (28) and (29) the acceleration of the center of 

gravity of the model and balance system will be obtained. The 

numerical solution shows that the maximum acceleration occurs when 

the natural frequency has the maximum value in its range, which is 

4 X 70 = 280 Hz. Then the maximum acceleration for y = 1. rad. 

■'o 

2 

is about 300,000 rad/sec . Piezoelectric accelerometers which can 
be sensitive to 10000 g can be used to detect the shuttle accelerations 
of this magnitude. 

In order to find the energy of the system, the sum of potential 
and kinetic energy should be obtained. Thus for a conservative 
system 

E = P . E + K • E (30) 

or 

E,. = %(K^.cp^ + (31) 

Eb = 

Therefore, using the equations (24), (25), (26) and (27) and 
substituting into the equations (31) and (32) one can calculate the 
total energy of the system. 

Considering the design of the balance system, two cross bars 
are designed as the lateral and torsional springs. These two bars 



126 


are designed such that the natural frequency in both perpendicular 

directions (called the principal directions) have the same value. 

Since the moment of inertias in these two directions are about the 

2 

same (65.7 and 65.5 Ibm-ft ), the stiffnesses, and thus the diameters 
of the bars, are about the same. In the case of bending lift and 
drag were estimated in the report of May 3, 1971. The projectionof 
these two forces can always be found on the principal directions. 

Now, instead of changing the wind direction, the balance system 
will be related, and the whole system will start to vibrate in the 
direction of the resultant of the lift and drag. This angle is 
generally not the same as the wind angle and is called the vibration 
angle . 

A computer program is written and attached to this report which 
calculates the lateral and torsional displacement, velocity and 
acceleration. Equations (24) through (29) are used in this program 
to predict the motion of the model in the wind tunnel. Results 
are computed for every 7.5 degrees of the wind angle over two 
complete cycles. A subprogram (Mohr) is used in the program which, 
for an input of the principal stiffnesses, produces the stiffness 
at the vibration angle and it simply uses the Mohr Circle. and 

in the equations (24) through (29) are obtained by the Mohr 
Circle in the direction of vibration. 

For every angle of the wind the corresponding lift and drag 
are generally different and for each case the corresponding stresses 
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and strains are obtained in the stiff bars. These stresses and 
strains are obtained from the following relations and are constant 

through the bars. 


M d/2 

O’ “ T 
max I 


(33) 


a 

max 

e => — — 
max E 


(34) 


where 

rt = maximum stress 

max 

e = maximum strain 

max 

M = moment at joint point 

d = diameter 

A 

T S 3 s - a the moment of inertia 
64 

E = Young's modulus 
The condition is shown in Figure (2) 



The moment M is the effect of lift and drag in the corresponding 


directions . 







ACCELERATION x 10 (rad /sec 
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ACCELERATION FOR f = 4. (sec”’'), ou = 10 (rad/sec) 
AND WIND ANGLE 7.5° 


* Possible Test Procedure By Dimensional Analysis : 

Considering the number of variables involved in the calculations, 
the question arises as to how many experiments should be performed 
to cover all possible variations of the variables and to indicate 
the critical point (or points)? Especially for a random type of 
input, such as wind velocity profiles the matter will be more 
complicated . 

Energy can be considered a characteristic which shows the 
states of the system. If energy picks up or goes to infinity at 
certain values of any of the variables, then the critical point is 
achieved. For a conservative system the energy can be obtained 
from the equations (31) and (32). Theoretically, energy is infinity 
if = u) or X|^ = U), that is if the natural frequency is the same 
as the frequency of vortex shedding. The conditions under which this 
may be achieved and the corresponding value of ou, should be inves- 
tigated experimentally. 

In order to lessen number of experiments for all possible 
variations of the variables, a dimensional analysis is performed as 
follows. The energy is generally a function of the following variables 

E = f(p, L, f^, f^, V, 0) (35) 

where 


* Book: Dimensional Analysis and Theory of Model, By Langhaar 

Wiley. 
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E = total energy = MLT 

_3 

p = density of air = ML 
L = characteristic length = L 

= frequency of vortex shedding = T ^ 

= frequency of lateral vibration = t” ^ 
f^ = frequency of torsional vibration = T ^ 

V = wind velocity = LT ^ 

9 = wind angle = 1 

For a given configuration 

= §(v> 0, p, L, shape of the model) (36) 

and since all these variables are listed in equation (35) f^ may 

be omitted. However, if the test may be continued for different 

models, then f^ should be listed as it is in equation (35). 

Regarding equation (35), one can varify that in order to have 

all the conditions and values of E, the variables f , , f , v and 0 

b t 

should vary within their range which requires a great number of 
experiments. By a dimensional analysis there is the possibility 
that a smaller number of parameters may be defined containing f^^, 
f^, V and 9 which can be used to determine all the necessary 
information. 

Eliminating f^ from the equation (35) and by a dimensional 
analysis the equation becomes 
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E 2 5 ^ 2 . 2_3v / 

- = f(L 0, L , V L ) (37) 

In the above equation f, and f are both considered to aflect the 

D t 

total energy. This is true in the case of a coupled vibration 

between the lateral and torsional inodes. In this case some 

vectorial form of f, and f^ will cause the critical condition. 

b t 

But if it has been proved that there is no coupling, then the 
equation (37) can be separated for lateral and torsional vibration, 
and in this case, the following equations will be used instead of 
equation (37) 




v\3) 


(38) 


E 


b 


P 


fbCL^fb^e. v\^) 


(39) 


However, referring to the more general case of equation (37) and 
calling 


E 

P 1 
L fb 0 = tT2 


t5. 2 

^ ^ = ^3 


2 3 

V L = TT/ 


(40) 


the equation (37) will become 
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TTi = f(n2» tt^) (41) 

Now, it can be seen that instead of varying f, , f , v and 0, one 

b t 

can vary tt 2 > and within their range and find all possible 
conditions . 

The test procedure should be as follows: 

1. Fix TT 2 and rr^ and find the variations of - TT 2 

2. Change one increment (depending on the accuracy of the data) 

and fix tt, and find the variation of rr, - tTo 
4 12 

3. Do this until the whole range of tt^ is covered 

4. Change one increment and repeat 1 through 3 

5. Do 4 until the whole range of is covered. 

If it is required that different models be tested in the wind tunnel, 
then the term f^ should be inserted in the equation (35) and this 
will cause another rr-term in equation (41) of the form 

TT. = L^f ^ (42) 

5 V 

In all of the above equations it should be noted that the 
characteristic length, L, and the density of air, p, are constant 
and are inserted to balance the equality. The characteristic length, 
L, can be any value such as the length or diameter of the booster 
or orbiter, but it should remain the same for all rr-terms. 
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Symbol List 

TL(I) = Lateral bending moment (due to drag) 

TT(I) = Tangential bending moment (due to lift) 

T(I) = The Resultant 

YZ = Initial rotation (.1 rad.) 

DT = Time Increment (varying; very small) 

TMAX = Max. time asked for (two complete periods) 

DTETA = Wind angle Increment (7.5°) 

XMMAX = Max. moment of Inertia (In the principal directions) 

XMMIN = Min. moment of Inertia (in the principal directions) 

XKMAX = Max. Stiffness (in the principal directions) 

XKMIN = Min. stiffness (in the principal directions) 

ALF = Angle of the resultant moment 
TETA = Wind angle 

BET = vibration angle = ALF + TETA 
XK = Stiffness in the direction of vibration 
XM = Inertia in the direction of vibration 
BAND2 = 

D 

D = Diameter of the stiff bars 

ST,SL,EPT,EPL = Stress Strain in the tangential & lateral bars 
E = Modulus of elasticity 

Y,YD,YDD = Displacement, velocity and acceleration of the center of 
gravity of the model and balance system. 
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When a fluid or gas flows past a cylindrical obstacle, the wake bcliind it is seen 
to contain numerous whirls or vortices airanged in a regular pattern, which was 
described and drawn accurately as early as the fifteenth century by Leonardo 
da Yinci. The whirls appear alternately on the one and on the other side of the 
cylinder, and they arc washed away in the wake when they have reached a certain 
eisc. At the moment that a vortex has reached its maximum siso, just before 
detaching itself from the cylinder, the velocity of the flow past that side of the 
cylinder is maximum and hence, by Bcnioulli's law, the pressure on that side is a 
minimum. The cj'linder thus experiences an alternating force in a direction 
perpendicular to that of the flow. 

In 1878 Strouhal published a formula for the frequency of this force or of this 
vortex shedding, based on observations only. This fomiula is now' usually repre- 
sented in a dimensionless form: S = JD/V, where S is the dimensionless “Strouhal 
number,” a constant of value 0.22 or about 1/4.5; / is the frequency, D the cylinder 
diameter, and V the stream velocity; the latter three quantities being expressed 
in any set of consistent units. This relation can be remembered b}' observing that 
one period of force oscillation (or the shedding of one pair of vortices) occurs during 
the time that it takes the stream to advance through 4.5 diameters. 

In 1911 von Karman made a stability analysis of the vortices and from it derived 
the geometrical pattern. This was the first theoretical investigation of the subject, 
and it is so important that the vortex wake now is generally referred to as the 
Ktirmdn wake. The intensity of the alternating force can be written, for dimen- 
sional reasons, as: .. . 



r.'-. yt 

■ Fk ~ Ckp — A 



.. where Fx is the Karman force; pV^/2 the stagnation pressure, -4 the projected 
area of the obstacle (perpendicular to the stream), and Ck a dimensionless coefficient, 
which may be called the Kdrman coefficient. It is expected that this coefficient is 
a function of the Reynolds Number. Calculations have been made by various 
authors, none verj' satisfactory, indicating that Ck =“ 1, in other words, that the 
maximum value of the Karman force is about the same as the stagnation force. 
Also, it is stated in the literature that for Reynolds Numbers over 3 X 10‘ for 
. cylinders, the Strouhal formula no longer holds, and the vortex wake is much 
. smaller, so that Ck ^ L The latter statement has been proved i:*'.'orrcct bv 
•recent observations of large industrial smokestacks vibrating violently at the 
correct Strouhal frcciucncy with Reynold.s Numbers as high as 7 X 10®. 

The abovo value of Ck is based on the a.ssumption that the cylinder or other 
obstacle docs not vibrate. If it docs vibrate; two effects have been noticed: First 
if the frequency of the cylinder vibration is in the neighborhood of the natural 


Strouhal frecpicncy (between SO and 120 per cent of it) then the cylijider motion 
forces the vortex shedding to take place at the cylinder frequency, anVl not at the 
Strouhal frequency, and second, the intensity of the force is greater for a cylinder 
vibrating at the Strouhal frequency than for a cylinder at rest. 

These two effects liavc been observed qualitatively only; no reliable numcricul 
experimental results arc available. The technical importance of airplanes is so 
very much greater than that of vibrating cylinders that not nnich attention h.as 
been p.aid to the intensity of the Karm.an force .so far. It is hoped that a compre- 
hensive experimental investigation of the Kdrmatr coefficient Ck as a function of 
(1) Re>molds Number, (2) amplitude ratio x/diametcr, and (3) frequency ratio 
Wr i/w-rou It "'ill carried out in the near future. 
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Tccluiicnl applications of the KdrmAn wnko have hecn teleiihonc lines, electric 
cross-country transmission lines, submarine periscopes, television broadcasting 
antennas, industrial smokestacks, singing ship propellers, and the failure of one 
large suspension bridge. In all c.7scs trouble is c.vpcricnced when the natural 
..frequency of the structure coincides more or less with the Strouhal frequency. 
A convenient measure for the relative intensity of the excitation in these various 
cases is the ratio between the exciting Karman force and the dead weight of the 
vibrating object: 

^ CK(PrVy^)Pl ^2 /A fV*\ 

T7 “ p,{r/i)DH r '"WW/ 

Here p//pt is the density ratio of the fluid and the vibrating Irody and the combina- 
tion V'/gD is the dimcnsionlcs.s Froude number. In general, a light, stiff structure 
of large surface area is most susceptible to those vibrations. Of the various ca.scs 
enumerated above, the first five have been described in the literature rcpcatedl}', 
and the la.st two are dealt uith in the papers referred to in footnotes 1 and 2. 
Here we discuss two recent cases of smokestack vibration. The first of these 
occurred in Monterrey, Calif., involving three stacks of 12-ft. diameter and 200-ft. 
height, made of welded steel. These stacks vibrated for several hours at the 
Strouhal frequency of about one cycle per second with largo amplitudes at a Rey- 
nolds Number of about 2.0 X 10*. Xo serious damage occurred, but to prevent 
recurrence, the stacks were provided with vibration dampers in the form of 
latticed stmts, placed between the three stacks and from them to the adjacent 
building. One end of each stmt was securely fastened, while at the other end a 
slipping clutch was provided. The slip force of this clutch was calculated to provide 
such damping of the motion that with an assumed Ck = 1 the rc.sulting amplitudes 
would be small. Since the installation of thc.se slip dampers two years ago, the 
vibration has not been observed. 

The next ca.se on record occurred in St. Clair near Detroit on several welded 
steel Stacies of 300-ft. height and IG-ft. diameter at a Strouhal frequency of about 
a cycle per second and a wind velocity of about 4.3 miles per hour; the Reynolds 
Number being about 7 X 10*. This resulted in a buckling of the stmcturcs in 
their upper part, .so that they had to be taken down to about half height and built 
up again. They are being provided with dampers in three guy cables located at 
120 deg. apart arouml the periphery near (he top of the stack. Where these cables 
come near the ground large flexible springs are put in them, which arc capable 
of changes in length of some two feet. The.se springs and the cable arc mounted 
in permanent tension, and a large truck-type automotive shock ab.sorbcr is in.stallcd 
across the springs. Again the damping of these absorbers is chosen so that their 
energy dissipation is somewhat greater than the expected energy input from the 
vortices, assuming a coefficient C*: = 1. 

In conclusion it can be stated that welded steel stacks of the size described or 
larger stacks arc un.safc, and that dampers or other means of suppre-ssing the 
vibrations are cs.scntial. 

• Kerr, Shannon, and Arnold, "The I’roblenw of the Singing Projicller,’* Inal. Mtch. Enjra. 
Proc., 144. 5 1 (in 10). 

• Farquhurson, •*Acio<lynumlr Stability of Su^iK'ni^ion Bridgi^s, the Tnrontii JCorrows Bridge,” 
{/niv. IPosA. Eng. Exp. SUi. Bull., No. JIO, I’urt 111. 




7.6. K^rm/In Vortices. \Yhen a fluid flows by a cylindrical obstacle, 
the w.'ikc behind the obstacle is no longer regular but in it will be found 
distinct, vortices of the pattern shown in Fig. 7.22. The vortices arc 



alternately clockwise and counterclockwise, arc shed from the cylinder 
in a perfectly regular manner, and arc associated with an alternating 
sidewise force. This phenomenon has been studied expcrimcnt.ally aud 
it has been found that there is a definite relation among the frequency/, 
the diameter of the cylinder D, and the velocity F of the stream, expressed 
by the formula 

= 0.22 (7.17) 

or the cj'linder moves forward by .about 4)-^ diameters during one period 
of the vibnation. It is seen that this fraction is dimcnsionlc.ss and that 
UuTcfore the value 0.22 is independent of the choice of units. It is known 
is the Slrouhal number. The eddy shedding on alternate sides of the 
♦.vlindcr causes a harmonically varying force on the cylinder in a direction 
Krpcndicular to that of the stream. The maximum intensity of this 
force can be written in the form usual for most aerodynamic forces (such 
4? lift and drag) as follows: . 

.Fk -- (Ck-VzpV*- A)' s\no>t (7.18). 

The subscript K is for K.4rmAn, Fk being the Kdrmdn force and Ck the 
Wiincnsionlc.ss) Kdrmdn-forcc coefficient. The value of Ck is not known 



141 


with grcftt accuracy; however, a satisfactory figure for it is Ck — 1 (good 
for a large range of Reynolds numbers from 10- to 10’), or in words; 
The intensity of the alternating force per unit of sidewise ])rojccted area 
is about equal to the stagnation pressure of the flow. 

The mechanism pf eddy shedding from a cjdindcr at rest is truly a 
self-excited one, because there is no alternating property in the approach- 
ing stream, and the eddies come off at the natural Stroulial frequency. 
Whereas for all previous eases discussed in this chapter the sclf-excitcd 
vibration in itself was dangerous, this is not the case here: we do not 
much care what happens in the stream; all that we care about is to know 
that a force of magnitude (7.18) is exerted on the cylinder. As a rule 
this is without much conscciucncc, and trouble occurs only if the sclf- 
cxcitcd frequency of eddy shedding (7.17) coincides with the natural fre- 
quency of the structure on which it acts. Then a resonance occurs which 
may be destructive. Objects on which this has been observed arc many 
and various: transmission lines, submarine periscopes, industrial smoke- 
stacks, a famous suspension bridge, a large oil storage tank, and tiny 
rain droplets. 

First consider a transmission lino of 1 in. diameter in a wind of 30 m.j).h. 
The Strouhal frequency (Ec;. (7.17)] gives 110 cycles per second. Vibra- 
tion of lines at such high frequencies and at small amplitudes have often 
been observed and have repeatedly ended in fatigue failures, llcsonancc 
obviously takes place at a high harmonic of the line in which the span i.s 
subdivided in many half sine waves, of the order of 20 to 30. Because 
of the high frequency and small amplitude of this motion, it has been 
found possible and practical to bring it under control by damped vibration 

absorbers. The simplest and most 
common construction, known o-s the 
“Stockbridge damper,” is shown in 
Fig. 7.23. A piece of steel cable with 
weights at its ends is clamped to the 
line. The cable acts as a spring and is 
roughly tuned to the frequency of the 
expected vibration. .\ny motion of 
the line at point A will cau.se relative 
motion in the pieces of cable, and the 
friction between its strands dissipates 
energy. The point A of attachment is so cho.scn along the line that it 
cannot coincide with a node of the motion, where the dam])cr would he 
useless. Since the length of the half sine w.ave varies between 8 ft. ami 
20 ft., a location at about G ft. from the point of support on a tower will 
fit most frequencies and wave lengths. This device, crude as it is, ha* 
proved to be entirely effective in protecting the lines against dam.agc 



Fio. 7.23. StocVbridgo damper for 
transmission lines, consisling of n piece 
of stranded steel cable about i:! in. long. 
cnrryinR c.sst iron sveifthts of about 1 lb. 
at caclt end. 


o 
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from the Kdrm^n-vorlcx vibration. An nttcinj)t to clc.sign one of (hc.se 
dampers for the gallo]>ing line of p.agc 209, where llic frccpienci' is 100 
times as slow’ and the arnpliluclc 100 times as large, would lead to weights 
of several tons, which is entirely impractical. 

Submarine periscopes, being cantilevers of some 20 ft. or more in length 
(wlien extended) and of diameters of the order of 8 in., have shown res- 
onance with Karinan vortices at .speeds of about 5 m.p.h. Since the fluid 
here is water with a high density p, the exciting force [Eq. (7.18)) is very 
large and the vibration is very severe, leading to a blurring of the view 
seen tlmough the device. A cure would consist in changing the cylin- 
drical cross section into a streamlined one, but the periscope tube must 
be free to rotate in its guide if views in all directions are wanted, and the 
constructional complications are so great that this has never yet been 
attempted. 

Steel industrial smokestacks regularly show Kdrm.4n-induced resonant 
vibrations at wind s))ccd.s of the order of 30 m.p.h. Brick or concrete 
slacks have not c.xhibitcd this phenomenon, and in steel stacks it has 
become worse when riveting was replaced by welding, thus decreasing 
the internal dampiiig of the stack. Recently (1953) a welded stack of 
IG ft. diameter and 300 ft. heiglit came to resonance at its natural fre- 
quency of about 1 cycle per second (in a wind of some 50 m.p.h.) 
with such violence that it buckled and developed a crack in tlie steel 
extending over 180 deg. of circum- 
ference. The damaged upper half had 
to be torn down and rebuilt, and the 
new stack W'as provided with vibra- 
tion dampers. Tliese were mounted 
in guy cables between (he top of the 
slack and the ground, .as shown in 
Fig. 7.24. There arc two large springs 

of several feet length, and the entire 
assembly of springs and guy wire 
is p\it in initial and jicrmanent ten- 
sion, so that the springs permanently 
are some 8 in. longer than when unstressed. Between the springs is a 
d.'ishpot in the form of a large automotive-truck shock absorber. When 
the stack starts vibrating, a relative motion occurs across the dashpot, 
tluis dis.sipating energy. The design is such that the amount of energy 
dis.sipated in the damper is about the same as that dissipated in the st.ack 
itself, thus doubling the natural damping, which is sufficient to prevent 
damage to the stack. 





Via. 7.24. Friction dninpcr built into .t 
r.uy wire of n sinokc.stnck. 



